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1. INTRODUCTION 

A number of algorithms are available in the literature that attempt to remove most of the effects of 
temperature from thermal multispectral data where the final goal is to extract emissivity differences. Early 
approaches include adjacent spectral band ratioing, broad band radiance normalization and the use of one 
band where emissivities are generally high (e.g., 11 to 12 pm) to determine the temperature (Salisbury, 
1992). More recent work (Salisbury, 1 992) has produced two techniques that use data averaging to extract 
temperature to leave a quantity related to emissivity changes. These two techniques (Thermal Log 
Residuals and Alpha Residuals) have been investigated and compared and appear to provide reasonable 
results. 


The analysis presented in this paper develops a thermal IR multispectral temperature/emissivity 
estimation procedure based on formal estimation theory, Gaussian statistics, and a stochastic radiance 
signal model including the effects of both temperature and emissivity. The importance of this work is that 
this is an optimal estimation procedure which will provide minimum variance estimates of temperature and 
emissivity changes directly. 

Section 2 discusses optimal linear spectral emissivity estimation and Section 3 is a summary. 


2. OPTIMAL LINEAR SPECTRAL EMISSIVITY ESTIMATION 


A stochastic model for the spectral radiance in the thermal IR using the Planck equation ignoring 
reflection and atmospheric contributions is; 
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c, First Radiation Constant 

c 2 Second Radiation Constant 

T Average background temperature in region 
e(A ) Average background emissivity in region 
X Wavelength 

Ar Change in background temperature (a zero mean random variable) 

Ae(A) Change in background emissivity (a zero mean random process) 

The term in braces in Equation 1 is the apparent emissivity, e^X). The model assumes samples 
are independent and identically distributed spatially. From Equation 1 it is clear that changes in 
temperature affect all wavelengths uniformly. Thus AT can be considered as a random process of rank 
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one with respect to wavelength. In a generalized sense it is also considered to be narrowband since all 
of the variation induced by temperature will fall along a single direction with a properly chosen expansion 
basis for the spectral radiance. The covariance of the apparent emissivity can be computed assuming that 
temperature and emissivity are independent as : 
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The second term on the right hand side of Equation 2 is (to within the temperature variance) known. 
Moreover in an absolute sense the second term more often is the dominant one (thermal IR phenomenology 
tells us this, especially in the daytime) since it is the one dependent on temperature. The first term on the 
right hand side of the equation is more difficult to specify. In order to describe it for a certain class of 
material an ensemble of measurements is needed for that class. Since only a very limited number of 
measurements exist for spectral emissivity in the thermal IR, a description of the first term for an arbitrary 
class is not available, in this development it will be assumed that the first term is stationary white. That 
is: 
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This assumption is very appealing for large regions because Targe regions usually contain many different 
classes of materials and it is reasonable to assume the second order statistics of the ensemble are white. 


2.1 Cancellation of Narrow Band Interference (Temperature) 

What follows here is a procedure for separating temperature and emissivity. What motivates the 
procedure is the observation that emissivity and temperature cannot be separated from a single 
multispectral measurement but a statistically based procedure can provide separation with acceptable error. 
Again writing Equation 1 but including sensor noise w(X): 

L(X,T) = jeU) + Ae(A ) + — j + < 4 > 


Since the means are assumed to be known or well estimated from the data, the mean radiance 
can be subtracted off along with a division by the blackbody function to obtain the following equation for 
observed variation in apparent emissivity: 

lt«ai = AeU) + Ci6 |_* — * v'U) (5) 


The quantity Ae(X) contains the relative emissivity behavior and it is this component that needs 
to be estimated. To estimate a realization of the random process Ae(X). the covariance of the apparent 
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emissivity is needed. The covariance of Equation 5 can be written as 
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Here N,/2 is the strength of the white sensor noise. It is recognized that the sensor noise term 
will not be stationary white due to the division by the Planck function. This will be ignored because modern 
sensors are generally clutter limited not noise limited, so a stationary assumption should introduce very little 
error. The second and third term of Equation 6 will be considered as known. Since it is assumed that 
the emissivity covariance is stationary white, the form of the apparent emissivity covariance is: 
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The estimation filter is required to be linear but not necessarily shift invariant with respect to 
wavelength. This implies that not only stationary but nonstationary random processes can be estimated 
as well. In fact the stochastic description of the temperature behavior clearly indicates that it is nonwhite 
and nonstationary due to its non Toeplitz form: 
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The optimal linear filter for estimating the temperature driven part of e* pp (A), h 0 (X,Ti), is obtained 
from the following integral equation (Van Trees, 1968,1971): 
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Equations of this form are fundamental to all linear signal processing. A series solution to this 
equation is obtained by using the eigenfunctions and eigenvalues of K AT (A,q). Since the temperature driven 
term is rank one there is only one eigenfunction and eigenvalue and the optimal filter can be written as 
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This filter, when applied to the stochastic portion of the signal, will produce an estimate of the 
correlated component driven by temperature. The filter projects the observables onto the temperature 
direction defined by €(A) / A. To produce an estimate of the stochastic portion of the emissivity, 
Ae (A) , the weighted projection is subtracted from the observations. The weighting is important because 
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if all the energy were to be taken out in the temperature direction the 'space* spanned by the emissivity 
would be dimensionally too small providing poorer estimates. To make the emissivity estimation equation 
simple in form the following are defined: 
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The optimal estimate of emissivity becomes 


A€(X) = Ae«P(A) - 


. |AeW(r|) ^(t,) dqj ^(X) 


(14) 


Equation 14 is the best linear minimum mean square estimate of the spectral emissivity which is the 
spectral apparent emissivity less the contribution from temperature variation. 


3. SUMMARY 

T emperature variations cause variations in the multispectral data that are highly correlated between 
spectral bands. Although many phenomenology based techniques have been developed to remove the 
variations due to temperature from the data, it has been shown in this paper that a formal development 
of the optimal estimator includes a temperature projection filter to remove the correlated variations due 
temperature leaving a direct estimate of emissivity changes. 
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